Spin currents and spin dynamics in time-dependent density-functional theory 
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We derive and analyse the equation of motion for the spin degrees of freedom within time-dependent 
spin-density- functional theory (TD-SDFT). Results are (i) a prescription for obtaining many-body 
corrections to the single-particle spin currents from the Kohn-Sham equation of TD-SDFT, (ii) the 
existence of an exchange-correlation (xc) torque within TD-SDFT, (iii) a prescription for calculating, 
from TD-SDFT, the torque exerted by spin currents on the spin magnetization, (iv) a novel exact 
constraint on approximate xc functionals, and (v) the discovery of serious deficiencies of popular 
approximations to TD-SDFT when applied to spin dynamics. 
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The dynamics of the spin degrees of freedom is respon- 
sible for such diverse phenomena as spin wave excitations, 
Bloch wall motion, spin-polarized currents, and spin in- 
jection and spin filtering; concepts and phenomena which 
are important, e.g., in the growing field of spintronics 
jij. The calculation of spin dynamics within density- 
functional theory (DFT) has consequently received much 
attention Q . The most popular DFT method for a first- 
principles treatment of the spin degrees of freedom is 
spin-density-functional theory (SDFT, see Refs. Q| and 
H for reviews.) SDFT has led to versatile and power- 
ful schemes for the calculation of, e.g., total energies, 
spin densities, and spin-resolved single-particle band- 
structures, but its traditional (i.e., ground state) formu- 
lation is applicable only to static situations. This situa- 
tion has changed with the advent of time-dependent DFT 
(TD-DFT) H , which has brought dynamical phenomena 
within reach of DFT. 

In the present paper we first derive the equation of 
motion for the spin magnetization from TD-SDFT, in- 
cluding exchange-correlation (xc) effects. We then show 
how this equation can be used to obtain information on 
the many-body spin current. Although one might think 
that a calculation of spin currents would require the more 
complex formalism of time-dependent cwrrent-density- 
functional theory (CDFT), we find that this is not en- 
tirely true: as a consequence of the continuity equation 
TD-SDFT suffices to calculate the spin currents in sev- 
eral cases of great of practical interest. 

Numerical applications of TD-SDFT, as of any other 
DFT, require knowledge of the xc potentials, which con- 
tain all many-body effects beyond the Hartree approx- 
imation. In traditional DFT many exact properties of 
these xc potentials are known, and greatly aid the con- 



struction of good approximations [|3|j4|], but the same is 
not true in the time-dependent case, where properties of 
the xc potentials are just beginning to be explored ||. As 
a byproduct, our analysis reveals a previously unknown 
exact property of the xc potentials of TD-SDFT, which 
strongly constrains suitable approximations. 

In TD-SDFT the fundamental variables are the time- 
dependent particle density, 
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n(r,t) = 53(*|*(r- ri )|*) 
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and the time-dependent magnetization (or spin) density 
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where ^ = ^(ri, .., r/y, t) is the many-body wave func- 
tion (spin coordinates are suppressed for brevity), /io = 
qh/(2mc) is the Bohr magneton, &i is the vector of Pauli 
matrices, and N the particle number. Here and below all 
operators are taken to be independent of time, i.e., the 
time dependence of expectation values results exclusively 
from that of the wave function. In TD-SDFT these ex- 
pectation values are not calculated with the many-body 
wave function but from the solutions of a noninter- 
acting Hamiltonian, containing suitably choosen effective 
electrostatic and magnetic fields v s {r,t) and B s (r,t), 
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Eq. (|3|) defines the Kohn-Sham (KS) Hamiltonian of 
TD-SDFT. The solution of [ihd/(dt) - H KS ]$ = 0, i.e., 
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the Slater determinant 5>(ri, .., rjy, t), reproduces the cor- 
rect particle and spin densities as functions of both time 
and position. By contrast, the current 
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j(r) = ^—Y. V^(r - n) + 5(r - ri )Vi, 
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is not among the basic variables of TD-SDFT, so that 
j(r,t) = (*|j(r)|*> ¥= <*|j(r)|*> = 3 KS (r,t). 

The equations of motion for the fundamental density 
variables describing the electronic degrees of freedom, 
n(r,t) and m(r,t), can be calculated directly from the 
commutator with the Kohn-Sham Hamiltonian. For the 
time evolution of m(r, t) one finds 



dm(r, t) 

dt 



V- J Ki (r,f) 
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m(r,t) xB,(r,t). (5) 



The second term on the left-hand side, which stems from 
the commutator with the kinetic energy, describes spin 
currents arising because the electrons carry their spin 
with them as they move around. Here 
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J KS (r,t) :=A*o> :<*|*i®ji(r)|3) 
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is the KS spin-current tensor, defined via the tensor prod- 
uct of the spin vector &i and the orbital current ji(r) 
of particle i 0. (We use a capital J to denote spin- 
current tensors, and a lower case j for the correspond- 
ing orbital currents.) Eq. (||) is a formally exact TD- 
SDFT representation of the time evolution of the spin 
degrees of freedom. Many-body effects enter Eq. (||) 
via the effective magnetic field B s , which is defined as 
B s (r, t) = B(r, t) + Ti xc (r,t), where B(r, t) is the exter- 
nal magnetic field, and B xc (r, t) the exchange-correlation 
magnetic field. 

Of course, the equation of motion for m(r, t) can also 
be derived from the many-body Hamiltonian. Since in 
the absence of relativistic effects the spin magnetization 
commutes with the particle-particle interaction, the re- 
sult takes the same form as in the Kohn-Sham system, 



dm(r, t) 
dt 



V-J(r,t) 
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m(r,f) xB(r,f), (7) 



up to the replacement of the effective by the external 
magnetic field, and the Kohn-Sham current 3 KS by the 
many-body current J. Eq. (fjj) is simply the continuity 
equation for the spin magnetization. Comparing Eqs. (||) 
and (Q) one finds that 

m(r, t) x B xc (r, t) = V • [J K5 (r, t) - J(r, t)\ (8) 



2mc 



for all times t and at every point r. This equation is 
the central result of the present analysis. We now pro- 
ceed to explore some of the consequences, practical and 
fundamental, of this result. 



(i) The right-hand side of Eq. (||) contains the differ- 
ence between the many-body spin current and its Kohn- 
Sham counterpart. This difference is in the following 
denoted 3 XC , since it is the exchange-corrrelation contri- 
bution to the full current. The longitudinal part of 3 XC , 
denoted J^ c , can be calculated directly from TD-SDFT, 
by integrating Eq. (||). It is interesting to note that this 
many-body correction to the Kohn-Sham current follows 
from SDFT alone, without having to use CDFT. 

J xc on its own is valuable information, but a full deter- 
mination of the many-body current requires knowledge of 
the transverse correction, too, and this can, at present, 
only be calculated within CDFT. For a large class of sys- 
tems, however, 3 XC is the most relevant many-body cor- 
rection. This includes all quasi one-dimensional systems 
(e.g., conducting polymers), in which only the longitu- 
dinal component of the spin current is of interest, but 
also three-dimensional systems in which the current flows 
entirely in one direction, e.g., perpendicular to an inter- 
face. A particularly important example is provided by 
magneto-resistive devices in which the spin current flow 
is perpendicular to the planes of a magnetic multilayer: 
such a current is purely longitudinal and can therefore 
be computed from Eq. (||). 

(ii) The left-hand side of Eq. (|J) is the torque locally 
exerted by the xc magnetic field on the spin configu- 
ration. The existence of this torque explains how TD- 
SDFT achieves a formally correct description of spin dy- 
namics via Eq. (|^) although the currents the spins are 
coupled to in that equation are the KS currents, and not 
the physical ones: the exchange-correlation torque (&) ac- 
counts for the difference between the many-body and the 
KS currents in Eq. (^J) . This torque arises from the com- 
ponent of B xc perpendicular to m(r, t), 
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(iii) Eq. (^) shows that for a spontaneously magnetic 
system (with B(r, t) = 0) the quantity —V • J(r, t) can 
be interpreted as the torque exerted by the current on the 
spin magnetization. There has recently been much in- 
terest in such torques in connection with experiments on 
current-driven magnetization reversals |J. Our Eq. (j|) 
provides a way to calculate these torques from the SDFT 
quantities H xc (r,t), m(r,i), and J KS (r,t). 

(iv) By integrating Eq. dq) over all space one obtains 



J d 3 rm(r,t) x B xc (r,t) 



0. 



(10) 



since the right-hand side vanishes by virtue of Gauss' 
theorem. Eq. ( |l0| ) has a simple physical interpretation: 
the self- consistent xc magnetic field cannot exert a net 
torque on the system as a whole. We refer to Eq. ( |To| ) 
in the following as the zero-torque theorem (ZTT). This 
terminology puts it on the same footing as the zero-force 



2 



theorem |^,|6| which states that the electrostatic xc po- 
tential does not exert a net force on the system. 

By integrating Eq. (||) over the volume of a finite sys- 
tem one finds instead of Eq. ( [To| ) 

/ Am(r,()xB IC (M) = f dS-3 xe (r,i), (11) 
J n Js 
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where f2 is the volume of the system bounded by the 
surface S. For a finite system the integral over m x B^c 
thus measures the net flux of xc spin current into and out 
of the system. 

(v) Finally, through Eq. (|l^) Eq. (||) provides a strin- 
gent constraint on approximations for B ac . Any such 
approximation which does not satisfy Eq. ( |l(i| ) is not con- 
sistent with the microscopic spin dynamics and the con- 
tinuity equation for the spin magnetization. In the pop- 
ular adiabatic local-density approximation (ALDA), for 
example, one takes as ones approximation the xc fields of 
static SDFT within the usual LDA, but evaluates these 
at the time dependent densities ||. Since in LDA ~B XC 
and m are by construction always parallel one has, in the 
ALDA, m(r,t) x B xc (r, i) = for all times t and at any 
point r, and Eq. ( |Io[) is trivially satisfied. The fact that 
the ALDA thus (fortuitiously) satisfies the ZTT may ex- 
plain the relative success LDA-based calculations have 
had in previous calculations of spin dynamics and other 
dynamical phenomena in TD-SDFT. 

On the other hand, from Eq. (||) one finds that within 
the ALDA J^ c (r, t) = 0. Since one does not expect that 
many-body corrections to the longitudinal spin current 
vanish generally, their absence in the ALDA must be 
considered a defect of this approximation. Furthermore, 
from the equation of motion (||) it follows that in the 
ALDA the xc magnetic field does not make a contribu- 
tion to spin dynamics, which is thus governed only by 
the noninteracting currents and the external magnetic 
field. In particular, for perfectly localized noninteract- 
ing moments (such that 3 KS = 0) and in the absence 
of an external magnetic field (B = 0), the ALDA yields 
dm/dt = 0, so that there is no spin dynamics at all. 
Within the linear-response approximation this problem 
of the ALDA has been noticed previously in Ref. 0. 
The absence of spin dynamics for localized moments is 
another deficiency of the ALDA. This conclusion is con- 
sistent with the original phenomenological derivation of 
the noninteracting equation of motion |ic(| , which shows 
that the spin dynamics of localized moments is driven by 
gradients of the magnetization, and these are not con- 
tained in the ALDA. 

Typical LDA-based calculations of spin-dynamics || 
avoid these problems by not proceeding entirely within 
the LDA (for example LDA is combined with constrain- 
ing fields which are not calculated selfconsistently, or it 
is used only to determine parameters in a model Hamil- 
tonian which itself is not of the LDA form). Any first- 



principles calculation fully within the LDA, however, is 
bound to run into the above problems. 

We thus now turn to a discussion of popular improve- 
ments upon the LDA. First, we note that any static func- 
tional that is invariant under the infinitesimal global spin 
rotation m(r) — » m(r) + 5(p x m(r) satisfies the ZTT, 
when used as input for an adiabatic approximation to 
TD-SDFT. An example is the exact-exchange (or OEP) 
method. Eq. (10) was derived only for exchange and cor- 
relation together, i.e., it is not guaranteed a priori that 
an exchange-only approximation satisfies it. However, 
since the exact exchange term is rotationally invariant, 
satisfaction of the ZTT is automatic. No similar result 
holds for a general time-dependent, e.g. retarded, xc 
functional. 

Next we consider gradient-dependent functionals. The 
static generalized gradient approximations of Ref. |Tl| l 
depend only on m z and Vm z . They also satisfy the ZTT 
identically, but at the price of making again B IC |j m. A 
gradient-dependent functional that depends on the full 
vector m and its derivatives was sketched in Ref. Q and 
another will be constructed below. The general form of 
such functionals is 

E xc = J d 3 re xc (n,\m\,Vn,Vm x ,Vm y ,'Vrn z ), (12) 

where the function e xc need not be specified explicitly 
for the present purpose. Using the chain rule for func- 
tional derivatives we write the derivative of the generic 
functional (Il2|) as 



<5m(r) 



5m k (r) 
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6m k [r) i(Vm,(x)) 



where the index k labels components of the vector m, 
and the unit vector in direction of component m k is 
denoted Ufc. The derivative in the first term on the 
right-hand side acts only on |m|, i.e., the terms Vm*, 
are held fixed while differentiating. By performing the 
derivatives of Vrrij with respect to mu and introduc- 
ing the abbreviations /(r) = <5/<5|m(r)| J d 3 r' e xc (r) and 
gfe(r) = <5/(5(Vmfc(r)) J d r'e xc (r) one obtains 



B xc (r) 



k 



u fc V • gfc(r), 



(14) 



where u m is the unit vector in direction of m. While 
the first term on the right-hand side is again parallel to 
m, the second, arising from the gradients in Eq. jl^), 
in general is not. Independently of the detailed form of 
the unspecified kernel e xc , we thus find that gradient- 
dependent functionals of all three components of m pro- 
duce a B IC that is not necessarily parallel to m. Such 
functionals thus have the potential to solve the problems 
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encountered by the LDA: when used as input for an adi- 
abatic approximation to TD-SDFT they can give rise to 
a nonzero B^ c (r, t) and hence, by virtue of Eqs. (||), (||), 
and (^|), to nontrivial spin dynamics for localized mo- 
ments, and to a nonzero J^ c (r,i). 

However, the ZTT constrains the dependence of the 
functional on the gradients of m„ m y and m z . In or- 
der to illustrate this with a simple example, relevant 
for the calculation of spin waves, we write B :EC (r) = 
/ dr' K [m, n] (r, r')m(r'), where the kernel K(r, r') is, for 
weakly inhomogeneous systems, a short-range function of 
r — r'|. We can then expand m(r') to second order in 
(r-r'), 

m(r') = m(r) + (r - r')V (g> m(r) + (r - r') 2 V 2 m(r), 

(15) 

and obtain 

B xc (r) = K 2 m(r) + K A V 2 m(r). (16) 

Here K.% and K4 are proportional to the second and 
fourth moments of the tensor K(\r — r'|), respectively. 
The ZTT is satisfied when these moments, which in gen- 
eral are tensors, reduce to scalars. The first term in 
( p"6| ) then recovers the Stoner approximation to the LDA, 
while the second is of the form of the terms appear- 
ing in the phcnomcnological Lifshitz equation of motion 
for slowly varying moments (cf. Eq. (6.98) of Ref. [tj3)- 
Work to construct explicit expressions for these moments 
is in progress. 

Up to this point we have considered conventional static 
SDFT only in so far as it provides the input for adiabatic 
approximations to TD-SDFT. We now take a brief look 
at static SDFT in its own right. The equation of motion 
for m in static SDFT is simply dm(r)/dt — 0. Hence 

V-J* S (r) = ^Un(r)xB s (r). (17) 

By comparing this with the corresponding many-body 
equation, which differs from it only through the replace- 
ment of B s by the external field B, and 3 K by J, one 
immediately finds 

-^-m(r) x B ac (r) = V • [J KS (r) - J(r)] . (18) 
zrac 

Since this is of the same structure as Eq. (j|) , our conclu- 
sions (i) to (v) hold in static SDFT, too, however now ap- 
plied to equilibrium currents and magnetizations, and not 
to their dynamical counterparts. Perhaps the most inter- 
esting result for static SDFT is the existence of a compo- 
nent of B IC perpendicular to m, since it implies that the 
prescription to take Ti xc locally parallel to m(r), which 
is employed in many recent calculations for noncollinear 
spin configurations |l3] is consistent with Eq. ( |l8| ) only 
to the extent that the right-hand side, i.e., the difference 
of the spin currents, can be neglected. 



In summary, we have derived and analysed the equa- 
tion of motion for the spin degrees of freedom within TD- 
SDFT. Our results have consequences for, e.g., the calcu- 
lation of spin currents in polymers and multilayers, the 
use of the ALDA in investigations of spin dynamics, pop- 
ular methods for treating noncollinear magnetism, and 
the construction of better density functionals. 
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